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Abstract 

We report on the measurement of the specific activity of ^^Ar in natural argon. 
The measurement was performed with a 2.3-liter two-phase (liquid and gas) ar- 
gon drift chamber. The detector was developed by the WARP Collaboration as a 
prototype detector for WIMP Dark Matter searches with argon as a target. The 
detector was operated for more than two years at Laboratori Nazionali del Gran 
Sasso, Italy, at a depth of 3,400 m w.e. The specific activity measured for '^^Ar is 
1.01±0.02(stat)±0.08(syst) Bq per kg of "^*Ar. 
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1 Introduction 



A 2.3-liter two-pliase (liquid and gas) argon drift chamber [T] was developed 
and built by the WARP collaboration as a prototype detector for WIMP Dark 
Matter searches with argon as a target [2] . The detector was operated for more 
than two years by the WARP collaboration at Laboratori Nazionali del Gran 
Sasso, Italy, at a depth of 3,400 m w.e. One important by-product of the op- 
eration of the prototype WARP detector was the precise determination of the 
^^Ar specific activity in natural argon. 

^^Ar and ^^Kr are two radioactive nuclides whose activity in the atmosphere is 
of the order of lOmBq/m^ and IBq/m^, respectively PH] . As a result of the 
liquid argon production process, they are both present in abundant quantities 
and are the two most significant radioactive contaminations in liquid argon. 
The two isotopes decay primarily by j3 emission, and their presence can limit 
the sensitivity of experiments looking for low energy rare events (WIMP Dark 
Matter interactions, neutrinoless double beta decay) using liquid argon either 
as a target or as a shielding material. 

^^Kr is not a pure (3 emitter, owing to the presence of a 0.43% branching 
ratio for decay with (3 emission on a metastable state of ^^Rb, which then 
decays by emitting a 7-ray of energy 514 keV, with a half-life of 1.01 /xs [S]. 
The coincidence between the /? and 7 emitted in a fraction of the ®^Kr decays 
may ease in some cases the task of determining experimentally the activity of 
^^Kr in low-background detectors. The determination of the specific activity 
of ^^Ar is intrinsically more challenging and not as widely discussed in the 
literature [511^ . A theoretical estimate is presented in Ref. [7j. 
In the last twenty years liquid argon technology has acquired great relevance 
for astroparticle physics applications. Several experimental techniques, em- 
ploying liquid argon as sensitive medium, have been proposed especially for 
rare events detection piSMTOllTTllTMMl] . For WIMP Dark Matter direct 
detection, the discrimination of nuclear recoils from the /3-7 induced back- 
ground plays a crucial role. The discrimination provided by the experimental 
technique must be sufficient to reduce the radioactive background below the 
very low interaction rates foreseeable for WIMP Dark Matter. A precise de- 
termination of the intrinsic specific activity of "^^Ar is therefore of significant 
interest for the design of WIMP Dark Matter detectors employing argon as a 
target. 



2 The 2.3-liter WARP detector 



The detector consists of a two-phase argon drift chamber with argon as a tar- 
get. Two-phase argon drift chamber was first introduced within the ICARUS 
program [5] in the framework of a wide-range study of the properties of noble 
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Fig. 1. Schematic view of the WARP 2 hter detector. 

gases. 

The drift chamber (see Figure [1]) is operated at the argon boihng point (86.7 K) 
at the atmospheric pressure of the Gran Sasso Laboratory (about 950 mbar) [T5] 
The cooUng is provided by a thermal bath of Uquid argon, contained in an 
open stainless steel dewar, in which the chamber is fully immersed. The pres- 
sure of the gas phase on the top of the chamber is naturally equalized to the 
surrounding atmospheric pressure. 

Ionizing events inside the liquid argon volume produce scintillation VUV light 
mainly at 128 nm. The scintillation light is shifted to wavelengths in the blue 
by an organic wavelength shifter (TetraPhenilButadiene) covering the walls, 
and collected by photomultipler tubes (PMTs) located in the gas phase and 
facing the liquid volume below. The 2-inch PMTs are manufactured by Elec- 
tron Tubes Ltd (model D757UFLA) and have a special photocathode that 
ensures functionality down to liquid argon temperatures. The quantum effi- 
ciency at the emission wavelength of TPB is about 18%. The material of the 
PMTs has been selected for high radiopurity: according to the supplier's spec- 
ifications, the total 7 activity above 100 keV is 0.2 Bq/PMT, dominated by 
the ^^^Th and ^^^U chains. 

A series of field-shaping rings surrounding the liquid phase superimpose an 
electric field of IkV/cm. The electrons are drifted toward the anode (located 
atop the chamber) and then extracted from the liquid to the gaseous phase 
by a local extraction field provided by a couple of grids. The electrons are lin- 
early multiplied in the gas phase by a second, stronger local field. The PMTs 
detect the primary scintillation light (directly produced by the ionizing event) 
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and also the secondary scintillation light (produced by the electron multipli- 
cation process in the gas phase). The PMT signals are summed and sent to 
a multi-channel analyzer recording the pulse height spectrum. The liquid ar- 
gon contained in the chamber is Argon 6.0 supplied by Rivoira S.p.A. The 
liquid argon is successively purified from electronegative impurities down to 
an equivalent contamination of less than 0.1 ppb of O2 by using the chemical 
filter Hopkalit™ from Air Liquid. The purity from electronegative elements 
is actively maintained by means of continuous argon recirculation through the 
chemical filter. 

The experimental set-up is located in the Laboratori Nazionali del Gran Sasso 
underground laboratory (3,400 m w.e. of rock coverage). The fiux of cosmic 
ray muons is suppressed by a factor of 10^ with respect to the surface (resid- 
ual fiux l.l/i/(m^- h), average muon energy 320 GeV [16]). The detector is 
shielded by 10 cm of lead, to reduce the external 7 background. 
The sensitive volume of the detector has the shape of frustum of cone and it 
is delimited by a stainless- steel cathode. The sensitive volume for the config- 
uration under analysis is 1.86±0. 07 liter. The density of liquid argon in the 
operating conditions (950mbar and 86. 7K) is 1.399 g/cm^ [15]. The sensitive 
volume is viewed by seven PMTs, whose responses have been equalized in 
gain. Daily calibrations ensure the long-term stability and the linearity of the 
response. The sensitive volume and the argon thermal bath are contained in 
a stainless steel dewar, 50 cm internal diameter and 200 cm internal height. 



3 Data analysis 

For the measurements described in this work, the electric fields were switched 
off and the chamber was operated as a pure scintillation detector. The gain 
of the PMTs has been set up to optimize the data acquisition in the typical 
energy range of the environmental 7-ray background, namely up to 3 MeV. 
The energy threshold for data acquisition is about 40 keV. The threshold used 
for analysis is 100 keV, in order to exclude events from electronic noise. The 
response of the detector to 7 radiation was studied using different 7-ray sources 
(^^Co, ^°Co, ^^^Cs) placed outside the chamber. The spectra obtained with the 
^^Co and ^^^Cs sources are shown in Figure [21 Typical values of the resolution 
observed with the calibration sources are a{E)/E = 13% at 122 keV (^'^Co) 
and a{E)/E = 6% at 662 keV (^^''Cs). The correlation between energy and 
primary scintillation light detected was linear within the range tested with our 
sources. The energy resolution of the detector can be described empirically by 
the following parametrization: 




(1) 
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Fig. 2. Energy spectra taken with external 7-ray sources, superimposed with the 
corresponding Monte Carlo simulations, (a) ^^Co source {E = 122 keV, B.R. 85.6%, 
and 136 keV, B.R. 10.7%), (b) ^^^Cs source {E = 662 keV). 

where ao = 9.5 keV, ai = 1.2 keV and a2 = 0.04. The three terms take 
into account the effects from non-uniform hght collection (02 term), statistical 
fluctuations in the light production (ai term) and electronic noise (ao term). 

The /?-7 spectrum in the detector (35 hours of live time) is displayed in 
Figure [31 The total counting rate is about 6 Hz (4.2 Hz above the analysis 
threshold). A simulation based on the Geant4 toolkit has been developed 
to reproduce and understand the observed features. The simulation takes into 
account the following two components of background: 



• External radiation from radioactive contaminants in the materials surround- 
ing the liquid argon volume (stainless steel, thermal bath, PMTs). Only 7- 
emitters are taken into account, that originate an effective 7-ray flux through 
the surface of the sensitive volume. This includes: ^ss-q g^j^^ daughters (es- 
pecially ^^^Rn dissolved in the external liquid argon), ^^^Th and daughters, 
^°Co and ^°K. The presence of these sources has been confirmed with a 
portable Nal 7-spectrometer inserted into the empty dewar. 

• Bulk contaminations in the liquid argon of the chamber. In this case, both 
/3 and 7 emitters are relevant. The most important contributions are from 
222Rn, 39 Ar and ^SKrE] 

The signal from internal ^^^Rn and its daughters can be monitored by counting 
the a decays of the isotopes ^^^Rn, ^^^Po and ^^^Po in the energy region 5—7 

^ Other radioactive isotopes of Ar, as ^'i'Ar and ^^Ar, are short-lived {T1/2 are 
35 days and 109 min, respectively) and their cosmogenic production is negligible 
in the underground laboratory. The long-lived ^^Ar (T1/2 = 32.9 y) is expected 
to be present in natural argon because of thermonuclear tests in the atmosphere; 
however, its concentration in '^'^''Ar is negligible, < 6 • 10~^^ g/g (at 90% CL) |18j . 
corresponding to less than 85 /xBq/liter in liquid argon. 
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Table 1 



Characteristics of the /3 decays of ^^Ar and ^^Kr. 



Isotope 


Half-life 


13 end-point 


(5 mean energy 






(keV) 


(keV) 


39 Ar 


269 y 


565 


220 


85 Kr 


10.8 y 


687 


251 



MeV. After a new filling with freshly produced liquid argon the chamber shows 
a total ^^^Rn decay rate of the order of 1—2 Hz. Due to the ^^^Rn half-life of 
3.8 days, the observed rate decreases to few tens of events/day four weeks after 
the filling. Therefore, j3 decays from the Rn daughters ^^^Pb and ^^^Bi can be 
neglected, provided the measurement is performed a few weeks after the filling 
of the chamber. The counting rate due to the decay of ^^C (dissolved in the 
liquid argon or located in the surrounding plastics) is estimated to be much 
less than 50 mHz. Most of the ^^C events occur close to the chamber walls. 
The main characteristics of the ^^Ar and ^^Kr (3 decays are summarized in 
Table [H Since both decays are classified as forbidden unique j3 transitions 
(AJ^'^ = 2~), the 13 spectrum is not described by the usual Fermi function. 
For the present work, we assumed the (3 spectra from Ref. [5]. 
The GEANT4-based simulation is used to generate normalized spectra Si{E) 
from the different radioisotopes, taking into account the energy resolution of 
the detector. The isotopes from the natural radioactive chains are treated 
independently. A fit of the experimental spectrum F{E) is then performed 
in the energy range from 100 keV to 3 MeV with a linear combination of the 
single components, i.e. 

F{E) = Y.w,-sm. (2) 



The coefficients Wi are treated as free parameters and represent the counting 
rates induced by the single sources. In Figure [3] we show the experimental 
spectrum, superimposed with the output of the fit (i.e. Y.WiSi{E)). The fit 
is satisfactory in all the energy range considered. The signals from the most 
important external 7-rays radioactivity sources and from internal contamina- 
tions are shown in Figure HI as derived from the analysis of the experimental 
spectrum. 



4 Discussion 

Figure H] shows that the energy region 2—3 MeV is dominated by interactions 
of 7-rays from ^32'jh daughters, the region 1.5—2 MeV by 7-rays from ^^su 
daughters, and the region 0.5—1.5 MeV by 7-rays from ^°Co and ^°K. Below 
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Fig. 3. Black histogram: background energy spectrum observed running the 2 hters 
detector deep underground at LNGS inside a 10 cm thick Pb shielding. The super- 
imposed red histogram is the result of a fit with a Monte Carlo simulated signal 
(see text). 

0.5 MeV the main contribution comes from j3 decays from internal contamina- 
tions of ^^Ar and ®^Kr; the two isotopes account for 65% of the total counting 
rate between 100 and 500 keV. 

The cosmogenically-originated '^^Ar contamination of '^^*Ar in the troposphere 
was measured in Ref. [3] to be (7.9 ± 0.3) • 10~^^ g/g; the quoted error was 
statistical onl}0]. Since the liquid argon used for the experiment is produced 
from the atmospheric gas, a similar ssAr/^'^^Ar ratio is expected to be present 
in our sample. 

®^Kr is mainly produced as fission product of uranium and plutonium. Its 
abundance in the atmosphere is of the order of 1 Bq/m^, corresponding to 
about 4- 10^^^ g(^^Kr)/g('^'^*Ar) in air. Nevertheless, the distillation procedure 
for the production of liquid argon substantially reduces the ®^Kr fraction. The 
residual ®^Kr in liquid argon may vary in different batches of liquid. 
In order to better show the Ar and Kr signals. Figure [5] displays the spectrum 
obtained from the experimental data after subtracting the fitted contribu- 
tion from the other sources. The single ^^Ar and ®^Kr contributions can be 
disentangled from the different end-point energies, 565 keV and 687 keV re- 
spectively. Since ^^Kr and ^^Ar decays populate the same energy region of the 



3 The 39Ar/=^atAr ratio was measured for dating purposes. The knowledge of the 
absolute ^^Ar specific activity was hence not necessary. 
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Fig. 4. The energy spectrum observed with the 2 hter detector can be reproduced 
by: (a) an external component dominated by interactions of 7-rays coming from U, 
Th, ^''Co and ^''K radioactivity of the materials surrounding the liquid argon; (b) 
an internal component dominated by ^^Ar and ^^Kr /? contaminations inside the 
liquid argon. 

spectrum their estimates are anti-correlated, as displayed in Figure El 
The specific activity of ^^Ar in liquid argon resulting from the analysis is 
1.41±0.02(stat)±0.11(syst) Bq/liter (la) (see below for the discussion of sys- 
tematic uncertainties). The corresponding ss^r/'^'^^Ar mass ratio is (8.0±0.6) • 
10-16 g/g^ with errors summed in quadrature. The result is in excellent agree- 
ment with the atmospheric determination of Ref . [3] . 

From the fit it is also obtained that the ^^Kr activity in the sample under in- 
vestigation is (0.16±0.13) Bq/liter (Icr). In a previous measurement performed 
with a sample of Argon 99.999% from Air Liquid, it was found a ^^Kr activity 
three times larger. This indicates that a non-negligible ^^Kr contamination 
may be found in commercial liquid argon samples. In this case, an additional 
fractional distillation could be required to reduce the radioactive background 
for the WARP experiment. 

The systematic uncertainties are summarized in Table [2l The dominating item 
is related to the energy calibration of the detector response: since the discrim- 
ination between ^^Ar and ^^Kr is based upon the (3 end-points, the fit result 
is sensitive to the energy calibration and resolution. The uncertainty on the 
energy calibration in the range of interest was evaluated to be 2% (Icr) from 
the meaurements with the 7-ray sources; the corresponding ^^Ar systematic 
error is 6.5%. The second important contribution is related to the uncertainty 
in the active volume of the chamber. The filling level can be determined with 
accuracy of about 1 mm, and the diameter of the teflon container with the 



8 



„ 0.009 
> 

i:^ o.ooe 

iZ 
D 

3 0.007 
0.006 
0.005 
0.004 
003 
002 
0.001 



"100 200 300 400 500 SOD 700 

Energy (keV) 

Fig. 5. Residual energy spectrum observed in the WARP chamber after subtraction 
of the expected contribution from external 7-rays. The spectrum is well described 
by a combination (dashed line) of /? decays of ^^Ar and ^^Kr. 

Table 2 



Systematic uncertainties on the Ar specific activity. 



Item 


Relative error 


Absolute error (Bq/liter) 


Energy calibration 


±6.5% 


±0.092 


Energy resolution 


±1.3% 


±0.018 


Sensitive mass 


±3.8% 


±0.054 


Total 




± 0.11 



reflector fixed on it is known with a precision of about 2 mm. This corresponds 
to an uncertainty on the sensitive mass of 3.8%. 



5 Conclusions 

The best estimate of the '^^Ar specific activity in the liquid argon is (1.41 ± 
0.11) Bq/liter, or (1.01 ± 0.08) Bq/kg of natural Ar, or (8.0 ± 0.6) • 10"^^ 
g(^^Ar)/g('^'^*Ar). The value is consistent with the previous determination by 
H. Loosli [3]. The uncertainty in our measurement is mainly due to systemat- 
ics. 

The liquid argon sample under investigation shows a contamination of ®^Kr, 
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Fig. 6. Confidence regions {la, 2a and 3a) for the ^^Ar and ®^Kr specific activities. 
Systematic uncertainties have been added in quadrature. 

0.16±0.13 Bq/liter (la). 
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